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In science we often use established code that was developed many years ago. Frequently the 
documentation is lost and many of us know the feeling of putting a number in the wrong column 
and trying to understand why the program is not working. However, these old codes are usually very 
useful, well tested and established. Rather than writing new code from scratch, which must first 
be tested against the established standard, it is desirable to put a new front end on the standard, 
commonly referred to as 'wrapping'. This paper will discuss wrapping of FORTRAN code using a 
selection of other languages to provide simpler user interface. 



I. INTRODUCTION 

In the ideal world data processing and modeling tech- 
niques are thoroughly studied and understood at the be- 
ginning of a research program. Then someone sits down 
and writes easily maintainable and extendible code with 
lots of comments throughout. Finally someone writes a 
detailed user manual including examples and what meth- 
ods were used, then confirms that everything does indeed 
work on all of the supported platforms. However, this is 
almost never the case. 

Normally, one is trying to understand some science and 
has little time to do more than write ad hoc code that 
frequently has a confusing interface and only is known 
to work for the parameters of the original problem. Af- 
ter time this type of software is grouped together in an 
analysis or modeling package that is fully functional, but 
lacks a real interface. Because of a steep learning curve 
to use the software a new user must have lots of interac- 
tion with the author to start to apply it to the current 
situation. This is a waste of everyone's time. The solu- 
tions to this problem are to undergo extensive rewrites 
of the code to simplify its user interface, or to wrap the 
code with new software which deals with all of the id- 
iosyncrasies of the various input files. This paper will 
discuss the latter option. 



II. SIMPLE WRAPPING 

Several cases of wrapping fall into a simple class where 
the underlying FORTRAN runs either off of a input file 
or a few command line parameters. This section will 
demonstrate an easy way to call such code. 

We will start with the absolute simplest program to 
call: hello world. This will be done with a very sim- 
ple FORTRAN program that requires no input or user 
interaction. 

program hello 

c 

write(*,*) 'Hello world' 
end 

Once compiled this will have either the name hello in 



Linux or hello . exe in DOS. For the rest of the paper we 
will assume that Linux is being used. To confirm that the 
program works we simply try (using a GNU compiler): 

7. g77 -o hello hello, f 
7o hello 
Hello world 

7. 

Even though this example is very simple it represents the 
vast majority of analysis and modeling software. The 
wrapping program in these cases can write the input file 
and call the FORTRAN using parameters that have al- 
ready checked to avoid errors in execution. 

Most require input either from a simple ASCII file 
(which can be set up by the wrapping software also) or 
a few command line parameters and be called in exactly 
the same way as our simple FORTRAN program. The 
examples in this section do not go over writing input files 
or parsing output files because the methods for calling the 
FORTRAN program is unchanged. 



A. Hello Perl 

In Perl 1,2 , and many other languages, a simple hello 
world type of program can be executed in a few lines. In 
this case the preferred function is system and the pro- 
gram is as follows (assuming that Perl is in /usr/bin) 

# ! /usr /bin/perl 
$command=" hello" ; 
system ($command) ; 

This can be tested on the command line 

7, hello.pl 
Hello perl 

% 

The system function gives control of the current thread 
to the function call and takes it back once the program 
executed is done. The exit status of the program exe- 
cuted (multiplied by 256) is the return value of system. 
There is also a similar function, exec, which starts the 
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system call but never regains control. Because of the na- 
ture of the two commands, system is the generally pre- 
ferred method since it provides a means to step through 
a series of programs in an orderly fashion. 

B. Hello Java 

Java 3 ' 4 is a significantly more structured language than 
Perl. Because of this, a simple 'hello world' program 
takes several more lines, while the basics are the same. 

import j ava . io . * ; 
class hello{ 

public static void main(String args[]){ 
String command= " hello " ; 
try{ 

Runt ime . getRunt ime ( ) . exec (command) ; 
}catch( java. io . IOException e){ 
// let it drop on the floor 

> 

} 

} 

Compiling and running this program results in 

°/ javac hello. java 
°/« java hello 

7. 

There is no output from the executed code. 

The lack of output comes from the exec command 
which grabs the input, output, and error streams from 
the program executed. To have the result of our program 
shown on the screen we need to ask exec for the process' 
output stream so we can redirect it to the console. This 
can be done by 

Process process=Runtime . getRuntimeO . exec (command) ; 
InputStream in_stream=process .getlnputStreamO ; 
InputStreamReader in_reader= 

new InputStreamReader (in_stream) ; 
Buf f eredReader in=new Buf f eredReader (in_reader) ; 
while (true) { 

output=in.readLine() ; 

if ( output==null ) { 
break; 

}else{ 

System. out. println (output) ; 

} 

} 

After the redirection lines are added execution results 

in 

7. javac hello. java 
°/« java hello 
Hello java 
% 

Much like the Perl example, this can be repeated with 
several different programs to allow for stepping through 
an analysis process. 



III. INTERACTIVE WRAPPING 

While the examples above represent the majority of 
wrapping needed, there are other cases which are impor- 
tant to discuss. Some other tasks are to deal with an 
interactive program or call FORTRAN subroutines di- 
rectly. Here we will discuss the former since the latter 
frequently requires an interface layer, such as C/CH — h, 
which calls the FORTRAN subroutines. Wrapping an 
interactive program happens when the original author 
spent the time to make a front-end to a program but 
new users want to use it through a GUI or with other 
programs in a seamless manner. 



A. KUPLOT and Perl 

The example here is taken from PDFgetN 5 which 
places a simple graphical user interface on analysis code 
written in FORTRAN. While several programs are called 
similar to the example above, part of the interface is plot- 
ting the data in various stages of the analysis. KUPLOT 6 
is a plotting package written using the PGPLOT library 7 . 
The three important subroutines arc to open a connec- 
tion to the outside program, send commands to it, and 
close the connection. Of course, if interactivity is not 
needed then the methods in the previous section can be 
used for any software. First we show how to open a con- 
nection to KUPLOT 

sub kuplot_open { 

my $cmd=File: :Spec->canonpath("$exepath/kuplot") ; 
(return) unless (-x $cmd I I -x "$cmd.exe"); 
my $pid = open2(\*kin, \*kout, "$cmd") 

II die ('KUPLOT connection failed ..'); 
kout->autof lush() ; 
kin ->autof lush() ; 

The important part of this subroutine is the open2 
method. This opens reading and writing connections to 
the same external code. The notation \*kin denotes that 
a reference is being used to pass the handle kin. Here kin 
is the output stream KUPLOT and kout is input stream 
of KUPLOT. In other terms, our program writes to kout 
and reads from kin in order to hold a conversation with 
KUPLOT. 

Next we send information to KUPLOT 

sub kuplot_send { 
my ($cmd)=@_; 
print kout "$cmd\n"; 
print kout "echo done\n" ; 
while (<kin>){ 

(last) if (/echo done/); 

} 

} 

A series of commands is given to KUPLOT followed by 
echo done. This is done so the output of the commands 
are parsed until the echo done is returned, which lets 
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the calling program know that all of the other commands 
were completed. 

Finally we close the connection to KUPLOT 

sub kuplot_close { 

print kout "exit\n\n"; 
close (kin) ; 
close (kout) ; 

} 

This is by far the simplest of the three subroutines. It 
tells KUPLOT to execute the exit command, then closes 
the input and output streams. 

B. BLIND and Java 

The example in this section is a program called BLIND 
developed for the Single Crystal Diffractometer (SCD) at 
the Intense Pulsed Neutron Source (IPNS) at Argonnc 
National Laboratory. Blind determines the orientation 
matrix, UB, 8 and lattice parameters consistent with a 
given set of Bragg peaks. Here the purpose of wrapping 
was to dramatically shorten the development time of vi- 
sualization and analysis software by simply using existing 
code and interacting with its results. Similar to the KU- 
PLOT example there are three basic steps: starting a 
process, interacting with it, and ending the process. The 
first step is fairly simple given the earlier Java example 

public Process startProcess (String command, 

String cwd) throws IOException{ 
// must specify a command 
if (command==null) return null; 
// get the working directory 
if(cwd!=null && cwd.length()<=0) cwd=null; 
if (cwd!=null){ 
File dir=null ; 
if (isDirectory (cwd) ){ 
dir=new File (cwd); 
return Runtime .getRuntimeO 

.exec(command,null,dir) ; 

} 

}else{ 

// call the shorter version of the command 
return Runtime . getRuntimeO . exec (command) ; 

> 

} 

This method allows for starting a process and specifying 
its current working directory, cwd . Without the extra 
parameter Java assumes that all processes are executed 
from the directory where the Java session was started. 

The process of interaction is similar to the previous 
example where things are read from and written to the 
process. Again, this method of interaction assumes a 
fair amount of predictability in how the program works. 
Since BLIND does not return a generic done message for 
each interaction, the possible results for each step and 
what to do next need to be coded into the wrapper. 

InputStream in_stream=process .getlnputStreamO ; 
OutputStream out_stream=process . getOutputStreamO ; 



InputStreamReader in_reader= 

new InputStreamReader (in_stream) ; 
OutputStreamWriter out_write= 

new OutputStreamWriter (out_stream) ; 
Buf f eredReader in=new Buf f eredReader (in_reader) ; 
Buf f eredWriter out=new Buf f eredWriter (out_write) ; 

output=this .readline(in) ; 
while ( output==null I I 

output . indexOf (" Input reflection from")<0 ){ 
if( output!=null && output. length()>0){ 
System. out. println(output) ; 

> 

output=this .readline(in) ; 

> 

this . writeline (out , "y") ; 
System . out . print In ( output+ " y " ) ; 

As in the simple Java example, the process input and out- 
put must be redirected by the caller to make the program 
work as expected. Lines of information are read from the 
process until the proper string is seen which needs actions 
taken. Then the appropriate answer is given to the pro- 
cess (and the console) to make it look like the program is 
being automatically run by Java. The reason for writing 
out the communication to the command line is to aid in 
debugging and provide a means to assure the user that 
the program is being run as expected. 

The readline and writeline methods are wrap- 
ping common functions in a convenience method. Since 
writeline is straightforward it is omitted. 

public String readline (Buf f eredReader in) 

throws IOException{ 
char buff[]=new char [BUFF_SIZE] ; 
String result=null; 
int go_back=0; 
int size=0; 

while( ! in.readyO ){ 

// wait until it is ready 

> 

for( size=0 ; size<BUFF_SIZE ; size++ ){ 
if(! in.readyO ){ 

if (size==0) return null; 
}else{ 

buff [size] = (char) in.readO ; 
String charac= 

(new Character (buf f [size] ) ) .toStringO ; 
if (charac . equals ( "\n") ){ 

go_back++; 

break; 

}else if (buff [size] ==-l){ 
go_back++; 
break; 

} 

} 

> 

if (size<go_back) return ""; 
result=new String(buf f , 0, size) ; 
return result; 

} 



4 



The reason for not using Java's default readLine method 
is that it stops reading only at the end of line character. 
Since BLIND has all of its input at the end of a line 
rather than on a separate line, some output is never fully 
written which causes the default readLine method to 
wait forever. 

Finally the process is ended by a simple 
process .waitForO . The connections to the input 
and output streams of the process are closed and 
memory freed automatically by Java. 

IV. SUMMARY 

This paper demonstrated methods for wrapping FOR- 
TRAN code. The techniques shown did not use any spe- 
cial interaction between the wrapping code and FOR- 
TRAN program so the same methodology can be ap- 



plied to wrapping software written in any language. Most 
modern languages have a method similar to the system 
command of Perl and exec command of Java. In cases 
where this, in combination with input and output files, 
is enough interaction between programs wrapping is very 
straightforward. Also discussed was the case of deal- 
ing with interactive sessions. While this is, in theory, 
straightforward some of the idiosyncrasies were demon- 
strated. 
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